Quantum capacitive phase detector 
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We discuss how a single Cooper-pair transistor may be used to detect the superconducting phase 
difference by using the phase dependence of the input capacitance from gate to ground. The proposed 
device has a low power dissipation because its operation is in principle free from quasiparticle 
generation. According to the sensitivity estimates, the device may be used for efficient qubit readout 
in a galvanically isolated and symmetrized circuit. 
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The interest in Coulomb blockade and accompanying 
quantum effects due to superconductivity has triggered 
a wave of research on new physics and future applica- 
tions. One basic device 0| is the single-electron transis- 
tor (SET) which consists of two small tunnel junctions 
having a sum capacitance of Cs. If the single-electron 
charging energy Eq = e^/(2Cs) dominates over temper- 
ature, Ec ^ ksT, the SET works as the most sensitive 
known electrometer. 

In order to gain advantage of the inherently large band- 
width [RsetCt.)-^ ~ 10 GHz of the SET charge de- 
tector, two new technologies have been developed where 
the SET is read using an LC oscillator. The RF-SET 
(radio-frequency SET) is based on the gate dependence of 
the differential resistance of a sequential tunneling SET, 
a property which modulates the Q-value of the oscilla- 
tor 1^. No other device than the RF-SET has been 
able to track dynamic single-charge transport at MHz 
frequencies, which is relevant especially from the point of 
view of characterization and eventual single-shot readout 
of superconducting qubits 

Because of the limitations due to the dissipative na- 
ture of the RF-SET, a technique called L-SET (Induc- 
tive SET), having low dissipation, has been developed 
recently With zero DC-voltage, a superconduct- 

ing SET, henceforth called single Cooper-pair transistor 
(SCPT), behaves as an energy-storing reactive compo- 
nent because of the Josephson coupling {Ej/2 is defined 
as the single-junction Josephson energy). 

Since the first energy band Eq of the SCPT grows 
approximately quadratically as a function of the drain- 
source phase difference 0, the SCPT behaves as an 
inductor when looked at from the source or drain. 
The effective Josephson inductance of a SCPT, L^^ = 
(27r/<I>o)^ d"^ Eq / dcj)^ , has a strong dependence on the (re- 
duced) gate charge Ug — CgVg/ (2e) if Ej/Ec ^ 1- Here, 
'I'o = h/{2e) is the flux quantum. In the L-SET schemat- 
ics, a charge detector is built so that the resonance fre- 
quency of a system of a SCPT and an LC tank depends 
on Ug. So far, we consider the L-SET the most promising 
method of sensitive and fast electrometry. 

In addition to resistance or inductance, capacitance re- 
mains in the group of linear circuit elements. The first 
energy band of a SCPT grows approximately quadrati- 



cally also with respect to the second external parameter 
Hg. (A corresponding statement holds also for a single 
junction, which has been suggested in Ref. Q as a ca- 
pacitance tunable by injected charge.) Therefore, while 
observed from the gate electrode, the SCPT looks like 
an effective capacitance CqB to ground. The operational 
principle of our proposal is based on the dependence of 
Ccff on the source-drain phase difference (f), 



dqg{4>y 



(1) 



where qnis the charge on the gate capacitor, as illustrated 
in Fig. |l| (a), (b). Note that in general qg ^ CgVg. 

Once coupled to a tank circuit, modulation of Ccff(0) 
can be used to read a phase difference (henceforth, called 
simply phase) reactively in a reflection measurement (see 
Fig. Phase is defined as the time integral of volt- 
age, = (27r/$o) /o Vdt. In the superconducting case, (j) 
equals the order parameter phase. 

Unlike any previous considerations of single-electron or 
single Cooper-pair devices, implementation of the "quan- 
tum capacitive phase detector" device proposed in the 
present paper is a generally fast and sensitive phase de- 
tector. We call the device a CSET, to emphasize its 
somewhat dual operation with respect to the L-SET de- 
vice. 

Hamiltonian for a SCPT symmetric in its Josephson 
energies is 

= ^cp(-|^ - - Ejcos{cp/2) cos{e), (2) 

where a term CgVg / 2 and terms having </> have been 
ignored. Ecp = (2e)^/(2Cs) is the Cooper-pair charging 
energy and Ej/2 is the Josephson coupling energy of the 
individual junctions. The phases are defined with the 
help of voltages Vi at points 1 and 2 in Fig. |l| = 
(2^/$o) /o - V2)dt and = (27r/$o) Jp Vidt. Here, 
the former is the difference and the latter is the sum of 
the phases over the two Josephson junctions. 

The eigenvalues for this Hamiltonian are given by the 
solutions to the Mathieu equation when Ug ^ n/2, where 
n is integer ||, |^] : 



FIG. 1: a) CSET circuit, b) An equivalent circuit. 



i?,(0,ng) = {k + l-{k + l)[mod2] + 2n,(-l)^ , (3) 

I 



where A1^(r, q) is the characteristic value A for even 
Mathieu functions |^ with characteristic exponent r 
and parameter q. Unavoidable asymmetry of tun- 
nel junctions in a real device is easily incorporated 
into Eg. |3|. By substituting Ej cos{(j>/2) in Eq. || by 
EjV (1 + c?^ + (1 — (P) cos(0))/2, we get energies of a 
SCPT whose individual junctions have unequal Joseph- 
son energies Ej{l + d)/2 and Ej{l — d)/2, where d ^ is 
the asymmetry parameter. Junction capacitances, how- 
ever, can be arbitrarily distributed. Fig. || illustrates the 
two lowest eigenvalues Eq and Ei with respect to the 
control parameters (j) and Ug. 

In order to calculate the observable capacitance Cog 
from gate to ground when the system is in the lowest 
eigenstate Eq, we first calculate qg — Cg{Vg ~ V2), where 
the island voltage is V2 = l/Cg{dEo/dVg). Using Eq. (1) 
we have 

where Cq is the quantum capacitance Cq^ = 
{d^Et))/{{2efdnl) due to the SCPT band structure. In 
the following analysis, the constant term Cg is neglected, 
because it is small compared with the shunting capaci- 
tance Cq. 

In order to get maximum performance of the phase 
detector, we take the operation point of the device such 
that the transfer function dCcS / d(f>, which is the deriva- 
tive of the capacitance modulation curves in Fig. 0, is 




FIG. 2: The lowest (|0>) and the first excited (|1)) energy 
surface for a symmetric SGPT (d = 0) with Ej / Ecp~3. 
The black line denotes the GSET operation point (f> ^ n and 

Vg{t) = VgO COStJot. 



maximized. This happens at (/> rather close to tt. As seen 
in Fig. ^, the transfer function increases rapidly at large 
Ej/Ecp- The price to pay for a high gain then is a lim- 
ited dynamic range. The external gate drive is taken as 
^git) = VgoCOSLUot. 

Asymmetry in the SCPT weakens the modulation con- 
siderably at high Ej/Ecp ^ S as seen in Fig. ^ (dashed 
lines). This is because asymmetry removes the degen- 
eracy at Ug = ±0.5 and (p = ±7r and smooths out the 
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strongest modulation. 
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FIG. 3: Capacitances Ccff calculated using Eq, 
2 fF and Ecp = l K. Solid lines are for a symmetric SCPT 
{d — 0), and the dashed lines for a slightly asymmetric SCPT 
(d — 0.1). The numbers denote the Ej/Ecp ratio. The circle 
and the square correspond to the operation points marked in 
Fig. ^ for the charge-phase qubit readout. 

As seen in Fig. ^ (b), the proposed CSET circuit is a 
series resonator with inductance L and total capacitance 
C = Co + Ceff. Resonance frequency / of the configura- 
tion is then tunable by Ceg. A shift of / due to a change 
of capacitance AC is 



/ = /o + A/ (X 



AC 



VC'o + AC 2C 



3/2 ' 



(5) 



and we find A///o ~ — AC/2Co. AC can be written in 
terms of a phase change Acj) as AC = (9Ccff/9(/))A0. For 
simplicity, we choose A^ = 1 radian. 

We assume that the resonator itself does not dissipate 
power so that its internal quality factor Qi = oo. The 
loaded Q-factor — 1 /Qi -f l/Qe is then set by the ex- 
ternal impedance Zg = 50 f2, so that Q = Qe = ojL/ = 
y/LjC/Zo. This is a good approximation with low Q. It 
follows that the phase d = aig{Vout/Vin) of the volt- 
age wave amplitude reflected from the resonator, as il- 
lustrated in Fig. 0, changes by tt in the frequency range 
of ~ /o/Q. At the resonant frequency, differential change 
of phase is larger, and by direct calculation one gets the 
maximum modulation Ai?/A/ = AQ/ fo, and 



Ai? 



2QAC _ 2ACVI 



(6) 



We assume operation safely in the linear regime, such 
that the gate charge amplitude is at maximum of the 
order e/2, and thus the gate voltage amplitude is Vgo ~ 
e/2Cg. The average energy stored in the capacitors CcS 
and Co is of the order E ~ K2^(Ccff+Co)/2 = e^C/iScf). 



Because it takes Q cycles for the (loaded) resonator to 
dissipate most of its stored energy, the reflected power 
flow is 



P - Efo/Q 



8gc2 ■ 



(7) 



This is about 5 fW with typical parameters in an ex- 
periment. The voltage amplitude of the modulation is 
AV = A-dVc, where the carrier voltage amplitude is 
Vc = \/2ZqP. The modulation is conveniently trans- 
formed into power units to enable later comparison with 
noise power. This leads to the definition of information 
power: 

(ATO! ^ (A^ ^ 
2Zo 2Zo ^ ^ ^ ' 

By combining Eqs. (|[), (0) and (||), we find 



2Cl ■ 



(9) 



Eq. (^) may be written in terms of the transfer function 
dC^s/d<j) = Cld^Eo/{Ae^dnld(j)), Co, and Z^. We also 
approximate Ccff ^ Co. It follows 



P.(A0) 



1 



QAirZoe 



( d'Eo 
^ \dn'id4> 



Cg_ 
Co 



(10) 



Sensitivity increases dramatically by increasing the 
Ej/Ecp ratio as is illustrated in Fig. ^. Eqp > kpT 
must be satisfied in order to keep the system localized in 
the lowest state Eq. At a fixed ratio of Ej/Ecp, max- 
imization of Cg/Co, and minimization of Zq yields the 
best sensitivity. 

In order to estimate device performance we take the 
operation frequency fo to be ~ 1 GHz. The individual 
component values are chosen to be easily realizable by 
present fabrication technology. Due to stray capacitance 
from bonding pads etc., it is difficult to achieve in prac- 
tice Co lower than ~ 0.15 pF. To have a series resonance 
at the fo with Q-factor of ^--^20, this would be accompa- 
nied by an inductance L ^ 160 nH. The inductance may 
be easily realized by using a commercial surface mount 
coil. It is to be noted that the capacitance Ccff is ^ 1/75 
of the capacitance Co assuming a large C^ ~ 2 fF. The 
large gate capacitance is easily implemented with a sim- 
ilar overlap junction as the two SCPT tunnel junctions, 
but using a much longer oxidation to create a highly re- 
sistive junction 

For phase sensitivity estimates, we calculated AC /Co 
resulting from a phase change Acf) — 1 rad, at different 
ratios of Ej/Ecp, and assuming Ecp = 1 K. They are 
listed in Table |. The values are calculated at the optimal 
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TABLE L Linearized relative change of capacitance corre- 
sponding to a A<j!> — 1 rad, for symmetric SCPT (parameters 
as in the text). 



Ej/Ecp 
AC/Co 



0.1 0.3 1 3 10 

..5e-5 6.6e-4 3.6e-3 l.le-2 4.3e-2 



Vgo in order to distinguish the carrier signal from the 
pre-ampHfier noise power within the bandwidth. For this 
reason, we have omitted the analysis of these points. 

In order to demonstrate an important application of 
the CSET, we discuss a circuit, where the CSET works 
as a detector for the state of the charge-phase qubit (see 
Fig. I). 



operation point of which maximizes the transfer func- 
tion, and at Vg ^ 0. As seen in Fig. ^, the information 
power is ~ 10 fW with Ej/Ecp ~ 10. 

By using a cryogenic HEMT amplifier with a typical 
noise temperature of 3 K that corresponds to spectral 
density of sn ~ 4x10"^"^ W/Hz, we find the phase sen- 
sitivity S0 = SAr/Pi(lrad) of the order 45 /irad/\/Hz 
for symmetric SCPT. Combined with a low-noise SQUID 
amplifier that would be close to the resonator cir- 
cuit, the value Zq could be lowered to a value of 1 
or below ||Tl[| and the noise temperature could be down 
by a factor 30. This would result in a sensitivity of 1.3 
/urad/VHz. 

Asymmetry in the SCPT junctions weakens the num- 
bers at high Ej/Ecp, as portrayed by the dashed line in 
Fig. § With Ej/Ecp = 10 and d = 0.1, the mentioned 
sensitivities would go down by a factor of 5. However, at 
high Ej/Ecp the junctions are naturally of large area, 
and hence relatively easy to fabricate with similar sizes. 
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Ej/EcP 



FIG. 4: Information power (Eq. corresponding to a phase 
modulation A(/> = 1 rad. Solid line: symmetric SCPT (d = 0); 
dashed line: slightly asymmetric SCPT (d = 0.1). Other 
parameters are as given in the text. 

For the estimates discussed here, we have assumed op- 
eration of the device only at "safe" values of Vg. That 
is, variation of the capacitance Ceff is continuous over a 
reasonable range of values around the operation point. 
As evident in Fig. ^, close to the degeneracy points, the 
differential change of capacitance is substantially larger. 
One needs, however, a relatively strong AC gate drive 




detector 



FIG. 5: Schematic of the CSET coupled to a charge-phase 
qubit The gauge-invariant phase differences 7^ and their 
orientations are marked with arrows. All Josephson junctions 
are taken to have same Josephson coupling energy E,j/2. $ 
is the externally applied magnetic flux. 

The basic idea of this circuit is that since the qubit 
states |0) and |1) correspond to distinct phases 0, they 
result in a different capacitance of the CSET phase de- 
tector. 

Due to a supposed large capacitance, Cg > 20 fF, 
phases with respect to ground at points A and B are 
taken as classical variables. Charging energy of these 
leads is neglected for the same reason. The structure 
having this high capacitance is naturally fabricated on a 
ground plane substrate. For a standard insulator thick- 
ness of 300 nm on top of a conducting ground plane, the 
loop structure would measure only tens of /im in size. 

We assume that the circuit can be described as two 
connected SCPTs (detector and qubit) and that the 
gauge-invariant phase 7^ over each of them is found by 
solving the equations 



71 +72 = -27r$/<i>o(mod27r), 



and 



9gfc'i(7i,K,i) _ dEk2{l2,Vi 



g2) 



dj2 



(11) 



(12) 



The former equation assures the correct 27r-periodicity 
of the sum of the phases around the loop. The latter 
equation states that the currents flowing through both 
SCPTs have the same magnitude. The energies E are 
defined according to Eq. (||). The band index k' for the 
detector is always null and the qubit index k takes values 
or 1, corresponding to the states |0) and |1). For this 
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example we assume for simplicity that both SCPTs have 
Ej/Ecp = 1- 

As an example of the operation we suppose that af- 
ter manipulating the qubit, it is in a superposition which 
then collapses into either |0) or |1) when the detector is 
turned on. The detector, on the other hand, stays in the 
state |0). In order to start the measurement, the qubit 
gate voltage is turned to approximately CgVg2/2e = 0.37. 
The detector is maintained at CgVgi/2e = throughout 
the operations. Thereafter, an externally applied mag- 
netic flux is ramped adiabatically to $/$o ~ 1/2. 

Fig. H illustrates the dependence of 71 with respect to 
the qubit state. Depending on the state of the qubit, the 
phase over the detector circuit 71 will thus become ei- 
ther or ~ —2 rad. According to Fig. ^, the capacitance 
Ceff then differs by ~ 0.4 fF, depending on the qubit 
final state. According to Eq. (^, the difference corre- 
sponds to an information power of —132 dBm, by using 
easily achievable values Co=0.15 pF, (5=20, /o=l GHz 
and Cg=2 fF. Signal-to-noise-ratio 1 is then achieved at 
a bandwidth of 20 MHz with a feasible SQUID first stage 
amplifier having a noise temperature of 0.3 K. Measure- 
ment is thus clearly possible in sub-microsecond regime. 
A reasonable junction asymmetry, d = 0.1 has a negligi- 
ble influence on the curves in Fig. ^. 

The proposed qubit circuitry has several important ad- 
vantages. First, low power dissipation means low rate of 
quasiparticle generation, which is essential for low back- 
action and fast recovery from the measurement. Second, 
since the circuit is galvanically isolated, it is free from 
external quasiparticle injection. 

The third advantage is the symmetry of the schemat- 
ics. Although in some sense the circuit is strongly cou- 
pled to the external Zq via the detector gate, thermal 
noise of Zq acts only as a common-mode signal. This is 
equivalent to say that the real part of the impedance seen 
by the qubit Re(Z) is very small. Asymmetry weakens 
the situation, but with a realistic, random asymmetry of 
about 20 % in the component values, we calculated the 
following figures in the schematics of Fig. ^ Re(Z) ^ 1 
mfl both at low frequency (/ < 1 GHz) which is rel- 
evant for dephasing, and at the level-spacing frequency 
(10 — 50 GHz) which affects relaxation. Especially since 
the qubit operations are naturally performed at the sad- 
dle point (j) = 0, Ug = ±1/2 similarly as in the original 
charge-phase qubit, the system is extremely well decou- 
pled from the environment. Note that the gate and flux 
operation leads of the qubit may be weakly coupled so 
that they do not contribute noise. 

Coherence time is then presumably limited by inter- 
nal 1// noise, which should be similar to existing qubit 
realizations. Its effect weakens as Ecp grows. At a con- 
servative value Ecp — IK, we estimate a dephasing 
time |l^ of 1 — 2 /is, which is comparable to the original 
charge-phase qubit. 

In conclusion, we have proposed a technique to mea- 



sure the superconducting phase difference by monitor- 
ing the effective capacitance between the gate of a sin- 
gle Cooper-pair transistor and ground. As a practical 
example, the readout of a charge-phase qubit using the 
technique was discussed. 
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pean Union. 




FIG. 6: Dependence of the phase over the detector 71 on the 
applied magnetic flux The thick (blue) lines correspond 
to the solutions of Eqs. ([l^) and (p^), when the qubit is in 
the state |0). The thin (red) lines are for the qubit state 
respectively. 
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